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Abstract 


Empirical modeling of plasma pressure and magnetic field for the quiet time nightside 
magnetosphere is investigated. Two models are constructed for this study. One model, referred to 
here as T89R, is basically the magnetic field model of Tsyganenko [1989] but is modified by the 
addition of an inner eastward ring current at a radial distance of ~3 Re as suggested by 
observation. The other is a combination of the T89R model and the long version of the magnetic 
field model of Tsyganenko [1987] such that the former dominates the magnetic field in the inner 
magnetosphere while the latter prevails in the distant tail. The distribution of plasma pressure 
which is required to balance the magnetic force for each of these two field models is computed 
along the tail axis in the midnight meridian. The occurrence of pressure anisotropy in the inner 
magnetospheric region is also taken into account by determining an empirical fit to the observed 
plasma pressure anisotropy. This represents the first effort to obtain the plasma pressure 
distribution in force equilibrium with magnetic stresses from an empirical field model with the 
inclusion of pressure anisotropy. The inclusion of pressure anisotropy alters the plasma pressure 
by as much as a factor of -3 in the inner magnetosphere. The deduced plasma pressure profile 
along the tail axis is found to be in good agreement with the observed quiet time plasma pressure 
for geocentric distances between -2 and ~35 Re- 
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Although the magnetosphere has been studied and surveyed extensively for over three 
decades since the advent of spacecraft, there are still outstanding tasks pertaining to the quiet tune 
magnetosphere. One of these is a quiet time magnetospheric model with specification of the 
magnetic field configuration and the associated equilibrium plasma pressure distribution. In the 
quiet time Earth’s magnetosphere, the inertial force pdv/dt is usually insignificant and therefore 
force balance should exist between magnetospheric plasma pressure and electromagnetic forces. 
An explicit configuration which resembles the magnetosphere and satisfies force balance has yet to 
be devised. 

Equilibrium of the magnetosphere has been studied theoretically in the two limiting 
situations of the inner magnetosphere where the Earth’s dipole field dominates and of the 
magnetotail where the field arises mainly from the theta-shape tail current system. The inclusion of 
the transition region between these two simplified situations has proved difficult [ Voigt and Wolf, 
1988]. In the inner magnetosphere, Sckopke [1972] studied equilibrium under the assumption of 
isotropic plasma pressure in a perturbed dipole field. This type of investigation is recently 
extended by Zavriyev and Hasegawa [1989] and Cheng [1992] with the more realistic 
consideration of anisotropy in plasma pressure as observed [e. g., Liu and Hamilton, 1992]. The 
magnetotail equilibrium has been studied analytically [Harris, 1962; Kan, 1973] as well as 
examined by numerical simulation [Toichi, 1972; Birn, 1989; Hesse and Birn, 1992]. 
Unfortunately, no global three-dimensional equilibrium solutions are known. Present global MHD 
simulations do not stabilize enough to provide them. Two-dimensional numerical solution, when 
convection is included, appears to evolve in ways which may preclude a steady state equilibrium 
[Erickson and Wolf, 1980; Erickson, 1984; Hau et al., 1989], albeit this may be merely an artifact 
of the two-dimensional limitation [Kivelson and Spence, 1988]. 

An alternative approach to obtain equilibrium configurations is from data. For example, the 
magnetic field models of Tsyganenko [1987, 1989], which give averaged magnetic field values 
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satisfying the static Maxwell’s equations, have been extensively used in relating magnetospheric 
regions to the ionospheric level and vice versa [Elphinstone et al ., 1990; Stasiewicz , 1991; 
Elphinstone and Hearn , 1992; PulkJdnen et al ., 1992]. However, early models of this type were 
badly out of force balance [Walker and Soutkwood , 1982]. Therefore, particularly pertinent 
questions related to any empirical magnetic field model are the following. What is the plasma 
pressure distribution in force equilibrium with the empirical field model? How is this plasma 
pressure distribution compared with actual measurements? It is quite evident that the availability of 
plasma pressure distribution in force balance with the magnetic forces in an empirical field model 
would further enhance its utilization. 

Spence et al. [1987, 1989] and Kan et al. [1992] have inverted the earlier magnetic field 
models [Tsyganenko and Usmanov, 1982; TsyganenJco, 1987] to obtain the equilibrium plasma 
pressure along the tail axis at midnight by assuming isotropic or nearly isotropic plasma. The goal 
of this paper is to extend these previous studies in deriving empirical magnetic field models for the 
quiet time magnetosphere, with the associated plasma pressure in force equilibrium with the model 
magnetic field and in good agreement with observations. Two features arc incorporated, namely, 
the presence of an eastward inner ring current and the occurrence of anisotropic plasma pressure. 
These two features are found to be essential in arriving at a realistic plasma pressure profile. This 
study treats the simple situation in which the dipole tilt angle is zero, making the configuration 
symmetric and easy to handle without altering significantly the results. 

the Modified Magnetic field models 

The magnetic field models of Tsyganenko [1987, 1989], referred hereafter as T87L (the 
long version) and T89, have a number of features which lend themselves to easy applications. The 
magnetic field in these models is given in an analytical form, with most terms bearing physical 
insights on the various current systems contributing to the total magnetic field at a given location. 
The coefficients of these terms are obtained through a least square fit to actual measurements and 
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are derived separately for different levels of the Kp index to allow for the variability of the magnetic 
field with geomagnetic activity. Since we consider here the quiet time magnetosphere, the activity 
level chosen for this study corresponds to the K p = 0 level. Modified versions of these field 
models form the basis from which magnetic field models are constructed here. 

The ring current is dominated by the gradient of plasma pressure. At the inner edge of the 
ring current, the plasma pressure decreases rapidly inward, producing an eastward-flowing ring 
current in the inner magnetosphere [Lui and Hamilton , 1992]. This ring current is featured in 
neither the T87L nor T89 model. Although this eastward- flowing ring current may not alter the 
magnetic field to a great extent, as will be demonstrated later, its absence affects tremendously the 
plasma pressure distribution required to maintain equilibrium with the magnetic forces. Without 
such an eastward-flowing current, the plasma pressure in force equilibrium with the magnetic field 
will have to increase monotonically inward, as indicated by the earlier work of Spence et al. [1989] 
and Kan et al. [1992], and is contrary to the observed plasma pressure distribution. Even though 
the required volume current density for the eastward ring current is small, the j x B force is 
significant since the field strength is high (about inversely proportional to the third power of the 
geocentric distance) in the inner magnetosphere. 

To include the inner ring current in the model, we construct its vector potential similar to 
the procedure used in the T89 model to represent the outer ring current, i. e.. 


A (mc) = C IRC pS r R C . (1) 

where Cjrc is a coefficient controlling the intensity of the inner ring current and 

= sj (P’ + (a IRC + 4®c) 2 » @) 

£irc = \A 2 + Po + fac ^irc) 2 > (3) 

^irc = 0-5[i + x (x 2 + L rc 2 ) ]• W 
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The location of the point in space is given in cylindrical coordinates by (p, 0, z) and in Cartesian 
coordinates by ( x , y, z). The physical significance of the nonlinear parameters are the same as 
given for the outer ring current by TsyganenJco [1989]. The radial scale length oirc relates to the 
geocentric distance for the inner ring current. The half thickness of the inner ring current region is 
represented by D 0 . The variation of this thickness with local time, i.e., the difference between 
dayside and nightside, is controlled by the function /iirc with its scale length Lirg We find that 
the flow reversal of the ring current from westward to eastward inside L ~ 3.5 can be reproduced 
reasonably well by adopting the following values for the above parameters: Circ = 760 nT, 
a IRC =1.5 Re, O q = 1.8 Re, TtRC = -0.2723 (same as the value for the quiet time outer ring 
current), and Lirc =1*5 Re- To maintain the total magnetic moment of the ring currents to be the 
same as the original ring current in T89 model, we have modified the outer ring current strength by 

Crrc = Corc ~ C IRC x (aiRC/flORc ) 2 > (5) 

where the subscripts RRC and ORC denote, respectively, the revised and original outer ring 
current parameters in T89. It turns out that the modification is extremely slight, only -0.2%. This 
model with the inner eastward ring current and a correspondingly modified outer ring current 
strength is referred hereafter as T89R 

Since the T89 and T89R field models give rather weak magnetic fields in the magnetotail 
[Stern and TsyganenJco , 1992; Rostoker and Skone, 1993; Peredo et al., 1993], we have also 
explored the T87L field model. To take advantage of the fact that the T89R model is better than the 
T87L model in the inner magnetospheric region but the reverse is true for the distant tail 
[ Tsyganenko , 1989; Donovan et al., 1992], we combine the two field models together by 
constructing a hybrid field model (referred hereafter as TH) with a transition parameter e such that 


B z ,TH = £ Sz,87L + (1 * £) Az,89R . 

£ = 0.5 - 0.5 sin ((2(x -xl)/(xu -jcl) -l)7t/2) , 


( 6 ) 

( 7 ) 
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where B z is the z component of the magnetic field along the midnight meridian, and subscripts 87L 
and 89R denote parameters from the T87L and T89R models, respectively. The parameters xl and 
xu are the lower and upper bounds of x for the transition. The transition parameter e is 0 at x = xu 
and 1 at x = xl, with a smooth derivative de/dx at the end points. The boundary points adopted 
here are xu = -5 Re and xl = -25 Re- We note that the region of interest in this paper is the 
equatorial plane along the tail axis in which the magnetic field has only the z component (for a 
dipole tilt angle of 0*). Therefore, the divergence of the magnetic field along the tail axis reduces to 
dB z /dz. It is clear then that the linear combination of the magnetic fields from the T87L and T89R 
models as expressed in (6) is also divergence free along the tail axis. On the other hand, if this 
hybrid model is to be extended beyond the tail axis, there would be a small non-zero divergence if 
the other field components are not modified accordingly. 

Figure la shows the profiles of magnetic field values in logarithmic scale from the four 
magnetic field models, i. e., T87L, T89, T89R, and TH, at the equatorial region along the 
midnight meridian. It can be seen that the four magnetic field profiles are quite different between 
x ~ -8 Re and -30 Re. The highest magnetic field values come from T87L, followed by TH, 
T89, and T89R. Although the inclusion of the inner ring current weakens only slighdy the 
magnetic field, this effect is manifested quite dramatically in the extremely low field region of the 
T89 model where the field magnitude goes below 1 nT. 

The volume current densities from these four models at the equatorial region along the 
midnight meridian are examined in Figure lb. The T87L model gives typically the lowest current 
densities among these models. Tail ward of x » -10 Re, the T89 and T89R models give almost the 
same values while the TH model shows noticeably lower current densities between x * -10 and 
-30 Re- The lower values from the TH model, in comparison with the T89 and T89R models, arc 
expected from the influence of the T87L model. Earthward of x - -10 Re, while the current 
densities from the T89R and TH models are almost identical, that from the T89 and T87L models 
have different values. Earthward of x » -8 Re, both the T89 and T87L models show a persistent 
westward current In contrast the T89R and TH models show a reversal of westward ring current 
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to an eastward ring current at -3.5 Re. The presence of this current reversal is modeled to 
represent the reversal noted from the CCE observations [Liu and Hamilton , 1992]. In comparison, 
we find the observed current densities appear to be higher than the model values at the current 
density peak but smaller than the model values at further downstream distances. In addition, the 
observed current densities are more variable, reflecting the filamentary nature of current much like 
the magnetotail current [McComas et al. , 1986]. 

The observed anisotropy of plasma pressure taken from Lui and Hamilton [1992] is shown 
in Figure 2 together with an empirical fit. The anisotropy is seen to be large in the inner region and 
decreases rather systematically at further downstream distances. The empirical fit to the anisotropy 
is a fifth-order polynomial function of x~ l , i.e., 

P±/P\\ - 1 = ao + a\ or 1 + a 2 x~ 2 + 03 r -3 + 04 r -4 + 05 x -5 , (8) 

with coefficients a 0 = -0.410554, a\ = -9.94369, a 2 = -86.9877, a 3 = -504.066, 04 = -1 1 10.73, 
and as = -847.912. The range of validity for this fit is -2.5 Re > x > -15 Re such that the plasma 
pressure is isotropic at x =* -15 Re- That the pressure should be isotropic downstream of x - 
-15 Re is indicated by the studies of Stiles et al. [1978] and Baumjohann and Paschmann [1989]. 

Given a magnetic field model, it is straightforward to compute the magnetic force j x B and 
determine the pressure force required to balance it. The radial profile of perpendicular plasma 
pressure can be obtained by integrating 

^=[jxB + {/> i -/’,Xt>-V)b], , (9) 


where the subscript r indicates the radial component of the vector quantity, b is the unit magnetic 
field vector B, j is the volume current density, P± and P\\ are the plasma pressure components 
perpendicular and parallel to the observed magnetic field direction, respectively. 
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Figure 3a shows the result from this computation for the T89R model with the plasma 
pressure at x = -34 Re taken to be 0.074 nPa. This initial value of plasma pressure is obtained by 
computing the plasma pressure required at the neutral sheet (z = 0) to balance the “asymptotic” 
magnetic field value for the tail lobe region of the T89R model at that distance. The “asymptotic” 
magnetic field value is the x-component of the field at z = 20 Re, which in this case is -13.6 nT. It 
is worth noting that since the plasma pressure is shown in logarithmic scale, a larger initial value 
for the plasma pressure will not affect the plasma pressure in the inner magnetosphere 
significantly. The plasma pressure is assumed isotropic at x = -15 Re and further downstream. 
The computed profile of the perpendicular plasma pressure is given by the solid trace. For 
comparison purpose, the profile of plasma pressure assuming pressure isotropy at all distances is 
given by the dashed trace. It is quite apparent that the assumption of isotropic pressure results in 
an underestimate of the plasma pressure at the inner magnetosphere, by as much as a factor of -3. 
This can be readily understood by noting that the plasma pressure anisotropy inside x =» -15 Re for 
the quiet time magnetosphere is in the sense of P± > P\\. With this sense of pressure anisotropy, 
the second term on the right of Equation (9) adds to the first one, thus requiring a larger pressure 
gradient for force equilibrium, as pointed out earlier by Lid and Hamilton [1992]. 

The computed perpendicular plasma pressure profile is compared with observations from 
Explorer 45, ATS-5, ISEE-2, and CCE in Figure 3b. The ISEE-2 data are taken from Spence et 
al. [1989] for the K p = 1~ condition, the CCE data from Lid and Hamilton [1992], the ATS-5 data 
from DeForest and Mcllwain [1971], and the Explorer 45 data from Spence et al. [1989] using 
measurements of 1-872 keV protons presented by Smith and Hoffman [1973]. The range of 
pressure values from Explorer 45 data indicated by the error bars actually spans the observed 
pressure taken during geomagnetically quiet to disturbed conditions and thus the lower limit should 
be used for comparison with the model calculation. As can be seen from Figure 3b, the model 
pressure profile is in general above the observed pressure profile, especially in the outer region. 
This is understandable and is to be expected because the observed pressure is not necessarily 
obtained right at the equatorial plane and thus is expected to be under-representing the actual 
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observed pressure at the equatorial plane. Inside 4 Re, the pressure determined from Explorer 45 
is noticeably lower than that from CCE values and the model calculation. This difference may be 
due to the upper energy threshold being only 872 keV for Explorer 45 measurements whereas that 
of CCE extends to >4 MeV. It is possible that a significant contribution to plasma pressure inside 
of 4 Re comes from particles above the 872 keV level. We have evaluated this possibility using 
several midnight passes of CCE during quiet conditions. We found the ions with energy 
>872 keV contributing typically ~50 - 65 % to the plasma pressure at geocentric distances inside 
of 4 Re. Since the disparity in the determination of pressure between the Explorer 45 and CCE is 
about a factor of 6, we conclude that the different energy passbands can account for a significant 
portion, but not its entirety, of the observed pressure difference. 

Overall, we find that there is a very good agreement between the observed and computed 
profiles of plasma pressure based on the T89R model. Another notable feature of the model 
pressure profile is its relative constancy between the downstream distances of 16 and 25 Re- This 
indicates the small magnitude of the j x B force in this region. Re-examining Figure 1 reveals the 
main reason for the decrease in the j x B force in the T89R model being the small value of the 
magnetic field in that region (Fig. la) since the volume current density is still substantial at those 
distances (Fig. lb). The more abrupt change in the radial profile of plasma pressure occurs in the 
downstream distances between x ~ -12 Re and x = -15 Re, as noted earlier by Spence et al. 
[1989]. Therefore, the j x B force is relatively large at the transition region between the dipole-like 
and tail-like field configurations. It is interesting to note that this region may correspond to the 
hinging point of the tail where the tangential stress acted on the magnetotail by the solar wind at the 
tail magnetopause is balanced by the attractive force between the Earth’s dipole and the tail current 
system [Siscoe and Cummings , 1969]. 

A similar plasma pressure profile can be obtained from the TH model, as shown in Figures 
4a and 4b. The plasma pressure based on the formula from Spence and Kivelson [1993] 


P (nPa) = We-*- 5 ™ + 8.9[xf'' 53 


( 10 ) 
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is also shown for comparison; the downstream distance x is in Rg. The assumption of isotropic 
pressure gives an underestimate of plasma pressure by as much as a factor of ~3 in the inner 
magnetosphere as before. The most significant difference between the two anisotropic pressure 
profiles in Figures 3 and 4 lies in the downstream region of x =» -10 to -20 Re- Whereas the T89R 
profile shows a relatively constant pressure in this region, the other profile shows a continuous 
pressure increase towards the Earth. This continuous rise, however, only contributes rather 
insignificantly to the plasma pressure Earthward of x « -8 Re- The continuous change in plasma 
pressure reflects a more gradual transition between the dipole-like and tail-like field configurations 
as well as a more ^-dependent pressure for the near-Earth tail region in this model than that in the 
T89R model. The empirical formula by Spence and Kivelson [1993] also compares reasonably 
well with the model profile. It gives a good fit to the plasma pressure within the range of x =» -2.5 
to -15 Re, but overestimates it inside this region and underestimates it further tailward. 


SUMMARY AND DISCUSSION 


We have investigated two empirical models for the specification of magnetic field and the 
plasma pressure for the quiet time nightside magnetosphere. The derived magnetic field models are 
extensions of the previous magnetic field models of Tsyganenko [1987, 1989] by the addition of 
an inner eastward ring current Furthermore, in constructing the equilibrium plasma pressure 
along the tail axis in the midnight meridian, we have taken into account the presence of plasma 
pressure anisotropy in the inner magnetosphere earthward of x = -15 Re- An empirical fit to the 
observed pressure anisotropy in the inner magnetosphere (between 2.5 and 15 Re) is also 
provided. These two new features affect profoundly the deduced plasma pressure profile. The 
first feature leads to an earthward reduction of the plasma pressure in the inner-most part of the 
magnetosphere (geocentric distances of < -3 Re) as observed. The second feature leads to an 
appropriate evaluation of the plasma pressure gradient. Without the consideration of plasma 
pressure anisotropy, the deduced plasma pressure may be under-estimated by as much as a factor 
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of -3. The derived perpendicular plasma pressure is found to be in good agreement with the 
observed values from ISEE-2, CCE, ATS-5, and Explorer 45. Therefore, the T89R and TH 
models provide quite realistic representations along the tail axis in the midnight meridian of not 
only the quiet time magnetic field but also the associated equilibrium plasma pressure distributions 
needed to provide the force balance with the magnetic forces in these models. The range of validity 
for the T89R and TH models in representing the magnetic field, the volume current density, and 
the anisotropic plasma pressure considered here is -2.5 Rg > x > -35 Re- The present result also 
supports the earlier findings from Spence et al. [1989] that the gradient of plasma pressure shows a 
relatively large change in the transition region between dipole-like and tail-like field configurations. 

A region called the inner edge of the plasma sheet has been introduced in the early studies 
of energetic particle environment in the magnetosphere [Vasyliunas, 1968; Frank, 1971]. It is a 
region characterized by an exponential decrease of electron energy density with decreasing radial 
distance and is generally located at geocentric distances of ~6 to 8 Re in the nightside. It is 
implicitly assumed by a number of researchers that the inner edge of the plasma sheet is associated 
with an earthward decrease of the total plasma pressure and theories on the closure of the large 
scale region 2 field-aligned current system [Iijima and Potemra, 1978] have been built upon this 
assumption (see, e.g., the review by Mauk and Zanetti, 1987). As pointed out by Mauk and 
Zanetti [1987], the observed plasma pressure shows a persistent earthward increase rather than a 
decrease in the geocentric distances usually ascribed to be the inner edge of the plasma sheet. In 
terms of total plasma pressure, there is no edge at those distances. The plasma commonly 
identified as constituting the plasma sheet in the outer magnetosphere (tail ward of the geocentric 
distance of -10 Re) gradually becomes the hot ring current particle population. The absence of a 
large scale plasma pressure decrease in the so-called inner edge of the plasma sheet indicates the 
necessity to modify the conventional mechanism for closing the large scale region 2 field-aligned 
current. An alternative means for the region 2 field-aligned current closure in an earthward 
increasing plasma pressure condition has been proposed by Lui and Hamilton [1992]. They have 
suggested, based on the work of Sato and Iijima [1979], that the correct sense of field- aligned 
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current will be generated if the VPj. vector is inclined towards the midnight meridian more than the 
VS vector in both the premidnight and the postmidnight sectors since the sense of the field-aligned 
current is determined by die sign of the triple product of B, VP±, and VS. Figure 5 illustrates this 
configuration graphically. In other words, a plasma pressure in the midnight sector slightly 
enhanced over that of the adjacent local time sectors will satisfy the requirement. 

There arc at least two aspects of this study which can be improved on in the future. One is 
to extend the calculation to other local times besides midnight and the second one is to extend it in 
the z-dircction to cover the region away from the equatorial plane. The former extension requires a 
better assessment on the variation of pressure anisotropy as a function of local time which is not 
presently available. The latter extension requires integrating the variations of both the 
perpendicular and parallel components of the plasma pressure in a plane. If these extensions were 
to be incorporated in the TH model, one must also address the extension of the divergence-free 
property of the magnetic field along the tail axis to other regions in the magnetosphere. This can be 
accomplished by considering a linear combination of the vector potentials of the T87L and T89R 
models instead of their magnetic fields. These improvements will be the subject of a future 
investigation. 
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Figure la Magnetic field profiles along the midnight meridian at the equatorial plane from four 
magnetic field models. 

Figure lb Current density profiles along the midnight meridian at the equatorial plane from four 
magnetic field models. 

Figure 2 Observed pressure anisotropy profile of the nightside magnetosphere and an empirical fit. 

Figure 3a Equatorial profiles of the perpendicular plasma pressure along the midnight meridian 
inverted from the T89R magnetic field model assuming isotropic plasma pressure (dashed 
trace) and adopting the pressure anisotropy shown in Figure 2 (solid trace). 

Figure 3b Equatorial profile of the perpendicular plasma pressure along the midnight meridian 
inverted from the T89R magnetic field model in comparison with measured values of 
perpendicular plasma pressure from several satellites. 

Figure 4a Similar to Figure 3a except that the perpendicular plasma pressure is inverted from the 
TH magnetic field model instead. 

Figure 4b Similar to Figure 3b except that the perpendicular plasma pressure is inverted from the 
TH magnetic field model instead. 

Figure 5 A schematic diagram to illustrate the relative orientation of the vectors VS and VSjl to 
give the correct sense of field-aligned current system for region 2. 
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